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The simulation of open quantum dynamics has recently allowed the direct investigation of the features of 
system-environment interaction and of their consequences on the evolution of a quantum system. Such 
interaction threatens the quantum properties of the system, spoiling them and causing the phenomenon of 
decoherence. Sometimes however a coherent exchange of information takes place between system and 
environment, memory effects arise and the dynamics of the system becomes non-Markovian. Here we report 
the experimental realisation of a non-Markovian process where system and environment are coupled 
through a simulated transverse Ising model. By engineering the evolution in a photonic quantum simulator, 
we demonstrate the role played by system-environment correlations in the emergence of memory effects. 

The simulation of quantum processes is a key goal for the grand programme aiming at grounding quantum 
technologies as the way to explore complex phenomena that are inaccessible through standard, classical 
calculators^'^. Some interesting steps have been performed in this direction: simple condensed matter and 
chemical processes have been implemented on controllable quantum simulators''^. The relativistic motion and 
scattering of a particle in the presence of a linear potential has been demonstrated in a trapped-ion quantum 
simulator^"^ that opens up the possibility to the study of quantum field theories^'^. The quantum Ising modeP° has 
been experimentally analysed under the perspective of universal digital quantum simulators'^ and very recently 
scaled up to hundreds of particles'^. 

This scenario has recently been extended to open quantum evolutions'^ marking the possibility to investigate 
important features of the way a quantum system interacts with its environment, including the so-called sudden 
death of entanglement induced by a memoryless environment'^. 

This interaction can destroy the most genuine quantum properties of the system, or involve exchange of 
coherence between system and environment, giving rise to memory effects and thus making the dynamics 
"non-Markovian"'^. Characterizing non-Markovian evolutions is currently at the centre of extensive theoretical 
and experimental efforts'^"'^. Here we demonstrate experimentally the (non-)Markovianity of a process where 
system and environment are coupled through a simulated transverse Ising model'" (see Supplementary Material). 
By engineering the evolution in a fully controlled photonic quantum simulator, we assess and demonstrate the 
role that system-environment correlations have in the emergence of memory effects. 

The paradigmatic description of quantum open dynamics involves a physical system S evolving freely accord- 
ing to the Hamiltonian and embedded in an environment E (whose free dynamics is ruled by H^). System and 
environment interact via the Hamiltonian H^^ '^, which we assume to be time-independent for easiness of 
description. While the dynamics of the joint state psE is closed and governed by the unitary operator 
U^ = ^-i{H^+H^+H^^)t^ the state of S, which is typically the only object to be accessible directly, is given by the 
reduced density matrix ps(t) = T^i^EiPsEit)}- For factorized initial states PseW — PsW ® ^^(0) and moving to an 
interaction frame defined by the free Hamiltonian of the total system, such reduced evolution can be recast in the 
operator-sum picture Ps(t)= ^^Kfi{t)PsW^l^ where {K^} is the set of trace-preseving, non-unitary Kraus 
operators of S that are responsible for effects such as the loss of populations and coherence from the state of the 
system'^. 

The roots of such decoherence mechanism have long been studied, together with their intimate connection to 
the so-called measurement problem and the implications of the collapse of the wave function^""^^. Numerous 
experimentally oriented techniques have been proposed to counteract decoherence'^'^^'^''. Yet, a complementary 
viewpoint can be taken, where the possibility to engineering structured environments and tailored S-E couplings 
is seen as a resource to achieve longer coherence times of the system^^ prepare entangled states, perform quantum 
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Figure 1 | The Sagnac interferometer used to perform our simulation. The red hne represents the high photon, whose polarisation is immediately 
detected. The low photon, depicted by the yellow and blue paths, goes through the set of gates acting on its polarisation and momentum, implemented by 
BS, HWP, and glasses plates The optical axis of HWP* is kept free so as to implement the last step of the simulation. The thin glass plates 0i, 02 and 03 
placed in the interferometer allow to set the phase of the environment evolution at each step. The Hadamard gates are implemented by setting 0i = 02 = 
03 = O[cf.Eq.(l)]. 



computation, and realise quantum memories. This calls for the ex- 
ploitation of memory- effects typical of a non-Markovian dynamics 
as a useful tool for the processing of a quantum state. 

Unfortunately, a satisfactory understanding of non-Markovianity 
is yet to be reached, which motivates the recent and intense efforts 
performed towards the rigorous characterisation of non-Markovian 
evolutions, the formulation of criteria for the emergence of non- 
Markovian features, and the proposition of factual measures for the 
quantification of the degree of non-Markovianity of a process^^"^^. 
Some of them have been recently used in order to characterize, both 
theoretically and experimentally, the character of system- environment 
interactions for few- and many-body quantum systems ^^"^^ Here we 
use a photonic setup to simulate a system-environment coupling ruled 
by a transverse Ising spin model (see Supplementary Material). Such 
interaction gives rise to the non-Markovian evolution of S, as wit- 
nessed by the non-monotonic behavior of its entanglement with an 
ancilla A that is shielded from the environmental effects. Our goal is 
the experimental investigation of the fundamental connection between 
system- environment correlations and non-Markovianity. The simu- 
lator that we propose allows for the implementation of various free 
evolutions of S and £, as well as the adjustment of their mutual coup- 
ling, thus making possible the transition from deeply non-Markovian 
dynamics, all the way down to a fully forgetful regime. 

Results 

Experimental setup and theoretical model. The simulator is experi- 
mentally realised using different degrees of freedom of the photonic 
device shown in Fig. 1, which consists of the concatenation of two 
Mach-Zehnder interferometers and a single Sagnac loop. The in- 
formation carriers are two photons (referred hereafter as "high" 
and "low"): the system S is embodied by the polarisation of the 
low photon (which can be either horizontal |H) = |0) or vertical 
I V) = |1)), while the environment E is encoded in the longitudinal 
momentum degree of freedom (the path) of the same photon (which 
will be right |r) = |0) or left |/) = 1 1)). The ancilla A is embodied by 
the polarisation of the high photon. High and low photons are 
emitted by the source of polarization- entangled described in the 
Methods. 

The ancilla is the key tool for our goals. Indeed, to investigate the 
emergence of non-Markovianity in the evolution of S due to its 
interaction with £, we use the method proposed in Ref. 27: we focus 
on the modifications induced by the S-£ coupling on a prepared 
entangled state of S and A. If the S-A entanglement decays mono- 
tonically in time, the dynamics of S is fully Markovian. Differently, if 
for certain time- windows there is a kick-back from E that makes such 
entanglement increase, the dynamics is necessarily non-Markovian. 



In fact, if the local action of the environment is no longer represented 
by a continuous family of completely positive maps, the S-A entan- 
glement is no longer constrained to decrease monotonically. This is 
evidence of the flow-back of coherence on the system and results in 
an increment of the S-A entanglement. 

As in other digital quantum simulators, the dynamics is approxi- 
mated by a stroboscopic sequence of quantum gates. Conceptually, 
the simulation consists of the forward evolution of S over discrete 
time slices^° according to a Trotter- Suzuki decomposition^^ of the 
total time propagator (see Supplementary Material). This approach 
is known to be effective for quantum simulation^ and is implemented 
here by the sequence of operations shown in Fig. 2a). The free evolu- 
tion of the environment is accounted for by the Hadamard gate 
= [g^ + g^) /\/2 {(T^rn IS the m = X, y, z Pauh matrix of qubit 
j = S,E and 1iJ is the Hadamard gate of qubit j = S,E), while 
H^ccgI. The S-E interaction is engineered by implementing the 
controlled-rotation G^^ = \0) {01^01^ + ®n\(p), which 

rotates the system according to the general single- qubit operation 
lZ^{(p)= cos (p^ — I sin (pG^y depending on the state of the envir- 
onment. This class of conditional operations is obtained from 
Hamiltonian generators of the two -qubit transverse Ising form (see 
Supplementary Material), which motivates our choice and is thus the 
class of system- environment interactions that is simulated in this 
work. Experimentally, we fix (p = 7i/4 so as to implement a con- 



ditional-((7j'H^) gate. By using the identity gI = WgIW, it is 
straightforward to prove that this gate is locally equivalent (via g^) 
to the composition of a controlled Hadamard and a controlled anti-Z 
(i.e. a gate applying gI only when E is in |0)), namely 
e^^ = C?^fc^CZfc^ = |0)(0|£®(7f + |l)(l|£(x)'H^. The extra needed 
to make the two conditional gates equivalent can thus be absorbed in 
H^y which thus becomes U^. 

This set of gates is experimentally realised in our photonic simu- 
lator as sketched in Fig. 2b). The gate is implemented by means of 
a beam splitter (BS) in conditions of temporal and spatial indistin- 
guishability of the optical modes. This scheme allows to evolve the 
input modes |r) and |€) at each step as 



|r)- 



1 



(|r) + e*|£)), 10- ^(k> 



(1) 



where the phases 0/ are varied by rotating thin glass plates intercept- 
ing one of the optical modes entering the BS. In order to perform a 
more general rotation of E it will be necessary to unbalance the out- 
put modes (using intensity attenuators) making the probability of 
occurrence of |r) and |€) unequal. 

The controlled gate CHk% is realised by placing a half-wave plate 
(HWP) on one of the two output modes with optical axis at 22.5° 
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Figure 2 | Scheme of principle of the simulation, (a) A maximally entangled state of A and S is prepared. S and E then evolve through the application of 
free single-qubit gates {H^ and H^) and joint two-qubits {CHknCZkn) ones. Time increases from bottom to top of the figure. Here nh(£) stands for the 
polarization of the high (low) photon, while for the momentum of the low photon, (b) Each block of the scheme is technically realised by the 
concatenation of a BS and two HWPs oriented at different angles. 



with respect to the vertical direction. The temporal delay introduced 
by this waveplate is compensated by another HWP on the opposite 
output mode with optical axis set at 0° with respect to the vertical 
direction. This implements the CZ gate, as shown in Fig. 2b). In order 
to assess non-Markovianity, ancilla and system are prepared in the 
maximally entangled state {(j)^)^^ = {\HH) + | VV)) ~/ Vl. This two- 
qubit Bell state is engineered using the polarization entanglement 
source of Ref. 32 (see Methods). The environment is initialised in 

\x)l= (^(x\r) + Vl-oc^\i)^^ (aeRis set by the transmittivity of 5Si 

[cf. Fig. 1] and can be varied by placing an intensity attenuator on one 
of the two output modes), a state endowed with quantum coherence 
as it is key due to the conditional nature of the dynamics that we 
simulate. The stability, modular structure and long coherence time of 
our interferometric setting allows for the repeated iteration of each 
block of gates {n^(g)l^) . 

As A does not evolve with the environment, the polarisation state 
of the high photon (depicted in red in Fig. 1) is immediately detected 
and only the low photon (in yellow in Fig. 1) goes through each 
sequence of gates at the various steps. A standard optical setup for 
the performance of quantum state tomography (QST)^^ is used to 
reconstruct the state of the S-A system after the evolution. The radi- 
ation is collected by using an integrated system composed of GRaded 
INdex (GRIN) lens and single-mode fibre^^, and is then detected by 
single-photon counters. For each step of the simulated dynamics, we 
measure the state of the environment by projecting the polarisation 
qubits on the states \HH), |HV), | VH), | VV). Finally, we trace out the 
degrees of freedom of the S-A system by summing up the corres- 
ponding counts measured for every single projection needed for the 
implementation of single-qubit QST. The Pauli operators for E, gen- 
erated after the first passage through the BSi, are measured using BS2. 
The same procedure is followed to perform the QST of the envir- 
onmental state at each step. The state of the S-A system is recon- 
structed in a similar way, by summing up the counts collected after 
projecting E onto |r) and 

Non-Markovian simulation. Figure 3 summarises the results ob- 
tained by running through the evolution of the overall system. In 
order to quantify entanglement we use the entanglement of for- 
mation EOF(SA)^^ between S and A, which is operatively linked to 



the cost of engineering a given state by means of Bell-state resources 
(see the Methods section). We are also interested in the correlations 
shared by the environment with the rest of the system, hence we 
evaluate the von Neumann entropy of the environment, defined as 
S{E) = —Ty[pe log2(p£)], which under the assumption of pure total 
ASE state quantifies the entanglement in the partition AS against E. 
Figure 3a) shows the experimental entanglement of formation at 
each time step (black-square points) against the results of a theo- 
retical model (red line) that, including all the most relevant sources of 
imperfections, deviates from the ideal picture sketched above. First, 
the BSs are not entirely polarisation-insensitive: for BSi and H (V) 
polarisation the reflectivity over transmittivity ratio R/T= 42/58 (45/ 
55), while for BS2 is R/T = 45/55 (55/45). Second, although the 
desired input state |<^^ )^ is created with high fidelity ( ^ 93%, see 
Supplementary Material), the entangled-state source generates 
spurious \HV) and \ VH) components, accounting for about 5% of 
the total state, which reduce the initial system-environment entan- 
glement to about 0.8. The inclusion of such imperfections makes the 
agreement between theory and experimental data very good up to the 
fourth step of our simulation, showing at least one revival of the S-A 
entanglement and thus witnessing the non-Markovian nature of the 
evolution^^. The fifth experimental point is significantly far from the 
theoretical behaviour because of the not-ideal setting of the phases 
(/ = 1, 2, 3). In fact, we have verified that their values are not 
completely polarization independent. This determines a slight dif- 
ference for the four contributions \HH)sa^ \HV)sA^\yy)sA^\yH)sA 
entering the state. This imperfection affects the performance of 
each Hadamard gate and becomes significant especially for the last 
step, where the cumulative effect of three gates should be consi- 
dered. Furthermore, the fifth point is affected by an uncompensated 
phase shift induced by the multilayer mirrors used in our setting [cf 
Figure 1], as explained in the Supplementary Information. None- 
theless, this point still reveals successfully the occurrence of a se- 
cond entanglement revival, thus strengthening our conclusions. 

Figure 3b) compares the experimentally inferred EOF(SA) (black 
squares) to the von Neumann entropy of the environment S(E) 
(green circles) quantifying the correlations shared between E and 
S-A. These figures of merit appear to be perfectly anti-correlated, 
thus giving evidence of a trade-off between the amount of entangle- 
ment that S and A can share at the expenses of S-E correlations. This 
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Figure 3 | Results of the simulation of the (non-)Markovian regimes, (a) Theoretical prediction (red line) against experimentally inferred data (black 
squares) for EOF(SA) as against the steps of the non-Markovian simulation. Two revivals of EOF(SA) are clearly visible, (b) Comparison between the 
experimental evolution of EOF(SA) (black squares) and S(E) (green circles) representing the correlations between E and S-A. The two evolutions are 
clearly anti-correlated. In order to measure the PauU operators needed to reconstruct the state of E at the fifth step, a further BS will be necessary, 
(c) Theoretical behaviour (red line) and experimental data (black squares) for EOF(SA) in the Markovian simulation. 



strengthens the idea that correlations with the environment play a 
fundamental role in this process, a point that has been addressed in 
Ref. 36 where it is shown that the establishment of system- envir- 
onment correlations is a necessary conditions for the emergence of 
non-Markovianity. This result can be bridged with our analysis con- 
sidering that the mixed system- environment state of Ref. 36 can be 
purified by enlarging the Hilbert space of the system including an 
appropriate ancilla. The anti- correlation between S(E) and EOF(SA) 
suggest interesting connections with the monogamous nature of quan- 
tum entanglement in such tripartite system^^. We have illustrated such 
connections in the Supplementary Information, where a theoretical 
simulation for up to 11 steps is reported and a quantitative analysis 
fof the entanglement sharing features in our system is provided. 

To reinforce our point even further, we now explore the Markovian 
counterpart of our simulation by replacing the unitary evolution of E 
with an incoherent map that resets the environment into the very 
same state at each step of the evolution. As the key role in the S-E 
interaction is played by quantum coherence, we chose to re-set the 
environment into a completely mixed state. This is realised by spoiling 
the temporal indistinguishability of the optical modes entering the BS 
by mutually delaying the |r) and |€) components. Intuitively, by mak- 
ing the state of the environment rigid, we wash out any possibility for 
system- environment correlations, thus pushing the dynamics towards 
Markovianity. This intuition is fully confirmed by the experimental 
evidences collected via QST: the black squares in Figure 3c) show a 
monotonic (quasi-exponential) decay of EOF(SA) (matching our 
theoretical predictions, red line), which is in perfect agreement with 
the absence of S-E quantum correlations as signalled by the positivity 
of the partially transposed S-E state. 

Discussion 

By simulating a non- trivial two-qubit coupling model, we have 
demonstrated the non-Markovianity of the evolution induced on S 
by a dynamical environment and a system-environment interaction 
allowing for kick-back of coherence. By adopting a witness that 
makes use of the effects that non-Markovianity has on entanglement, 
we have explored experimentally the link between the emergence of 
non-Markovianity and system- environment correlations. The next 
step in this endeavour will be the experimental proof that non- 
Markovianity can be used as a resource for the advantageous proces- 
sing of information^^'^^, such as the preparation of interesting states^°, 
along the lines of previous studies on state engineering and informa- 
tion manipulation through Markovian processes^^"^^. Our setup will 
be particularly well suited for this task, in light of the effective control 



over both system and environment that can be engineered both in 
space (thanks to the modular nature of the setting) and in time. 

Methods 

Source of entangled states. The black box shown in Fig. 1 of the main text represents 
the source used for the generation of the hyperentangled state encoded in the 
polarization and path degrees of freedom of the photons produced by spontaneous 
parametric down-conversion. Here polarization entanglement, corresponding to the 

state ) = — /= (|HH) + 1 VV)), was realized by spatial and temporal superposition 
of photon pair emissions occurring with equal probability, back and forth, from a 
single barium borate (BBO) type-I crystal under double excitation of a vertically 
polarized UV CW pump beam. A suitable rotation of the polarization on one of the 
two possible emission directions of the photons was then applied. The source that we 
have used is such that the two photons, belonging to the degenerate BBO emission 
cone, are also path entangled, as explained in Ref 32. For our pourposes we need to 
select a pair of correlated directions, belonging to the emission cone surface, along 
which two photons travel, obtaining in this way the high and low photons shown in 
red and yellow in Fig. 1. 

Entanglement of formation. For two qubits, the entanglement of formation is given 
by'' EOF = h I 1 , with h{x) = -x log2 x - (l - x) log2(l - x) and C = 

max{0, — >^2 ~ -^3 ~ ^4}- Here, {/I,} is the set of eigenvalues (arranged in non- 
increasing order) of the Hermitian matrix R = ^~p^fp where 
~P= {(^y®(^y)p* {(yy(E)(^y) and p* IS the conjugate of the density matrix. 
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